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Abstract 

A generalization of the Gibbs-von Neumann relative entropy is proposed based 
on the quantum BBGKY [Bogolyubov-Born-Green-Kirkwood-Yvon] hierarchy as 
the nonequilibrium entropy for an N-body system. By using a generalization of 
the Liouville-von Neumann equation describing the evolution of a density super- 
operator, it is demonstrated that the entropy production for an isolated system is 
non-negative, which provides an arrow of time. Moreover, following the procedure 
of non-equilibrium thermodynamics a master matrix is introduced for which a mi- 
croscopic expression is obtained. Then, the quantum Boltzmann equation is derived 
in terms of a transition superoperator related to that master matrix. 

Key words: quantum statistical mechanics, nonequilibrium and irreversible 
thermodynamics, kinetic and transport theory of gases 



1 Introduction 



According to the mechanicistic interpretation of the physical world, the basic 
laws of nature are deterministic and time reversible. However, at the macro- 
scopic level we observe irreversible processes related to energy degradation 
which generate entropy. How do we reconcile the 'spontaneous production of 
entropy' with the time reversibility of the microscopic equations of motion?. 
This is a problem not solved in a clear and definitive way yet, thus it still 
constitutes an open problem. Attempts to address this have been made based 
on coarse-graining procedures [l]-[3]. In a previous work [4] we showed that in 
the framework of the classical BBGKY- hierarchy which explicitly incorporates 
the correlation between all the particle clusters in the system, macroscopic ir- 
reversibility arises without coarse-grainig procedures. All that is needed is a 
coherent definition of the entropy given through a generalization of the Gibbs 
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relative entropy in terms of the distribution vector in addition to the gen- 
erahzed Liouville equation. Based on the fact that the generahzed Liouville 
equation introduces a trace preserving transformation in the set of distribu- 
tion vectors, we obtained a non-negative entropy production for the isolated 
system which manifests the existence of a classical time arrow by appealing 
to the underlying microscopic dynamics only. 

Here we develop the quantum version of this previous analysis. The quantum 
counterpart of the distribution vector is the density superoperator evolving 
according to the generalized Liouville- von Neumann equation. Moreover, we 
propose a generalization of the Gibbs-von Neumann relative entropy in terms 
of the density superoperator, in the spirit of the generalized Gibbs relative 
entropy postulate. As in the classical approach, here, the generalized Liouville- 
von Neumann equation introduces a trace preserving transformation aswell. 
Prom this we obtain a non-negative entropy production for the isolated system, 
providing an arrow of time and giving a microscopic basis to the second law of 
Thermodynamics. After obtaining that, through a thermodynamic analysis, 
the quantum Boltzmann equation follows in a natural way. 

We have structured the paper as follows. In section 2, we introduce the Hamil- 
tonian dynamic of the N-body system obtaining the generalized Liouville-von 
Neumann equation. In section 3, we compute the entropy production. Section 
4 is devoted to the non-Equilibrium Thermodynamic analysis of the system 
obtaining the quantum Boltzmann equation. Finally in section 5, we empha- 
size our main conclusions. 



2 Hamiltonian Dynamics 

Let's consider a dynamical system of N identical particles whose Hamiltonian 
is given by the N-particle Hamiltonian 



where H{j) denotes the individual energy of the j-th particle and (p {j, k) the 
interaction energy between the j-th and k-ih. particles. The state of the system 

is completely specified at a given time by the N-particle density operator 
TV) which evolves according to the Liouville-von Neumann equation 



(1) 



[5] 




(2) 
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In the Liouville representation, where the density operator is represented as 
a vector in the Liouville space [7]-[10], the Liouville- von Neumann equation 

becomes 

where repeated Greek indices mean summation and the Liouvillian L is defined 
through 



= La,,sD^,T , (3) 



La(S,^5 - [H^s hi " -^7/3 ■ (4) 

However, a strictly equivalent alternative description of the state of the system 
can be given in terms of the set of n-particle reduced density operators [5] , [6] 

P={l>(°),L>W, ,L>(^)} , (5) 

where T> constitutes a density superopcrator whose components are the n- 
particle reduced density operators D^'^\l^ ....,n) defined through 

^'"'(l.--.'')-(^,„,S,^„,^'^'(l......^) (6) 

and D^'jj — Sa/s {i-e. D^^^ coincides with the unit matrix). These operators 
satisfy 

^j,Jn<-Hi,....,n)} . (7) 



The superopcrator V can be represented alternatively as a (A?" -|- 1) x (iV + 1) 
block diagonal matrix in an abstract space 



JV -| 

P^^-T^^L'H, (8) 



n=0 

where the (N -|- 1) x (A^ + 1) matrices In are characterized for having zero 
everywhere except for the n-th row and the n-th column and satisfying 

In -^p ^np (9) 

which defines the inner product in this abstract space. Moreover, in view of 
Eq. (7) 

N -1 

= Tr {d^"^} = 2-. (10) 

On the other hand, the dynamics of the D^"") follows after performing the 
operation Tr on both sides of Eq. (2). Thus, we obtain the so-called 

(n+l,...,iV) ^ ^ ^ 

quantum BBGKY-hierarchy of equations which can be written in a compact 
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way as a generalized Liouville equation 

ih^Vit) =CQV(t) . (11) 
at 

At this point, it is convenient to introduce the projection operators V and Q, 
its complement with respect to the identity (i.e., VC + QC — £), which give 
the diagonal and non-diagonal part of the generalized Liouvillian C, respec- 
tively. This projection operators are defined according to 

N N 

VC=Y.Y. ^Pn^n ®^n® (12) 
n=0 p=0 

where L*^"^ is the n-particle Liouvillian 

L^:L - {H^Spy - H^pas) (13) 
and H^^^ is the n-particle Hamiltonian. Thus, it becomes 

N 1 

P£ P = ^ -X„ ® L>(")) . (14) 



n=0 



In addition, when the inner product involves two densities operators U and 
V, it should be understood as 

A(DB = A^pBp^ . (15) 

On the other hand, 

N N 

Q£=^T„®^5p(„+i)IV {Tp®A(f)}, (16) 

n=0 p=0 ^ 

where we have defined 

Ail,, = i^aX. - Fi%s) (17) 



and 



Therefore, 



F(^)(l,....,p) = Xi0(j,p) . (18) 



Q/:0P = V-X„0 Tr (a("+i) Z}("+iH . (19) 

^oUl (n+l)^ J 



Hence, Eq. (11) can be rewritten as 



ih^Vit) -VCqV^QCqV . (20) 
at 
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According to Eqs. (12)-(18), one of the components of Eq. (20) will read 

= Tr {Fir^)^/^) - F^+^^i^i"/^)} . (21) 

Without the right hand side term in this equation, the evolution of D*^"^ would 
be unitary. In the next section we will show that irreversibility is manifested 
in the dynamics of the system when the adequate description, i. e. in terms of 
the set T> of all the n-particle reduced density operators is used. 

3 Non-equilibrium entropy and irreversibility 

As the expression for the non-equilibrium entropy we propose 

S = ksTr {V {\nVo- \nV)} + So 

iV 1 N 

^ksYl -^n O E ® (in - In D^P^) + So 

n=0 p=0 

^ksY. -^^"^ (in - In + So , (22) 

n=0 ^• 

a functional of D^'^\ generalizing the Gibbs-von Neumann entropy postulate 
[11], [12], based on the fact that the superoperator V determines the state 
of the system. This entropy is analogue to the relative or KuUback entropy 
[17], [18]. Here, ks is the Boltzmann constant, Sg the equilibrium entropy and 
Vo is assumed to be the equilibrium value of V which corresponds with So, 
satisfying CVo = 0. It is known that S* is a concave functional of V, and will 
never decrease under a linear stochastic transformation like the one induced 
by Eq. (11) [14]-[18] 

V{t + t)^{U + tC) © V{t) , (23) 
with U being the unit superoperator 

N 

U = EXn^Xn^C/^'^l") (24) 

n=0 

and 

^(n|n)^ E Io.p,,s\a){P\®\^){5\ , (25) 

a/3, 7(5 

where /a/3,75 = 5^p5aS is the unit tetradic [19]. Therefore, the entropy we 
propose is coherent with the second law according to which S increases in 
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irreversible processes such as the relaxation to equilibrium from an initially 
non-equilibrium state. To prove this we compute the rate of change of S which 
can be obtained by differentiating Eq. (22) with the help of Eq. (11) 



where 
and 



= (7i + (72 > , 



ai^ -ih'^kBTr{VjCQV(\nVo-lnV)} , 



(26) 



(27) 



a2 = -ih-^kBTT{QCQV{\nVo-\nV)} . (28) 
Eq. (26) constitutes the entropy production corresponding to the relaxation 
passing from a non-equilibrium state to equilibrium. By using Eqs. (12)-(19) 
in addition to the cyclic invariance property of the trace, one has 



1 

ai = -ih-^ks E -T^^"^ ^^"^ (lnL»^") - lnL>(")) 

n=0 ^• 

N 1 

= -ih-^kB J2 -1 [^^"^ ^^"^] © (ini^i") - lnL>(")) 

n=0 ^• 

N -, 

n=0 ^• 

iV 1 

— 1 r),' L J 



(29) 



and 



a, = -i/l-iA;^ E ^ Tr {a("+i) © L'^^+^U © (lnZ}(") - lnL>(")) 

n=0 
N 1 

^-irtsE-T '^^ {[F^"+'\/^^"+'^]}0(lnD(")-lnL>(")) 

^(inD^'^)- In 



1 

ih-^kB E -r Tr If 



(n+l) 



(30) 



These are zero at equihbrium when = D^^\ but in any other case these 
are not necessarily zero because neither the density commutes with 

and Dl") nor D^"^ commutes with H^^\ in general. Hence, according 
to our previous discussion related to the concave character of the entropy 5* 
and its consequences, a should be non negative. This provides a time arrow 
inherent to the irreversible macroscopic manifestation of the dynamics of the 
system. 

As known due to the unitary character of the time-evolution operator related 
to the Liouville-von Neumann equation giving the dynamics of the N-particle 
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density operator, the von Neumann entropy is a constant of motion. Never- 
theless, as we have seen here, this is not the case for the generahzed Gibbs-von 
Neumann entropy, based on the fact that the dynamics of the density super- 
operator given through the genarahzed Liouville-von Neumann equation is 
not unitary. In addition, as compared to other definitions of entropy, Eq. (22) 
does not contain any coarse-grainig. In the hght of these facts not contained in 
other definitions of the entropy in the literature the strength of our definition 
(22) resides. 



4 Non-equilibrium thermodynamics analysis 



Here we proceed as in Ref. [20], which constitutes the standard way of Non- 
equilibrium thermodynamics. Thus, for small deviations from equilibrium we 
assume that V and Vg are diagonal in a common basis. According to this, the 
entropy (22) becomes up to second order 

S = —kB^^ {5VV-^5V} + So , (31) 
where 5T> — T> — Vg. Therefore, the entropy production is 

where the thermodynamic force 

X = SS/SV = ~kB (p-^SV+SVV-^) (33) 

is conjugated to the flux dV/dt. Prom Eq. (32) the phenomenological relation 

&D 

— ^TZQX (34) 

is inferred, where 

N N 

7^= ^$^2-„(8)2-p®i?("lf) (35) 

is a phenomenological supermatrix which due to the Hcrmitian character of 
the force X , we choose Hermitian on the right-hand side entry. This implies 
that 

D(n|p) _ D*("b) (nf'\ 



7 



where the superscript * stands for the complex conjugate. Using this property, 
one can replace X in Eq. (34) by 



N 

X^-kBV;^5V =Y.^^n®X^'^\ (37) 



n=0 

a nonHermitian form. Furthermore, due to the fact that T) is normalized, from 
Eq. (34) we infer that Tr(7^0 A") = 0, which leads to the following constraints 

E^S = 0. (38) 

a 

By substituting the phenomenological relation (34) into Eq. (32) we reach the 
expression 

AT -1 N N 

ih-' ^ Tr {[x("),L>("+i)] f("+i)} = ^ ^X(")0i?('^IP)0X(P) , (39) 

a kind of Green-Kubo relation, also obtained by using Eqs. (29) and (30) in 
the near-equilibrium approximation. 

Additionally, TZ should be Hermitian as predicted by the Onsager symmetry 
relations, which reflects in 

p("-lp) _ p*(pl") (An\ 

This leads to the constraint 

E<i = 0, (41) 

a 

and to the hermiticity in the left-hand side entry 

p(w|p) _ p*(n\p) fAr,\ 

With Eq. (37), Eq. (34) may be written in the form 

dV 

— = -Mq5V, (43) 

where the master matrix is defined by 

M = - kenV-^ . (44) 
According to the Onsager relations, it follows that the master matrix satisfies 
MVo = -kBll = -ke'R) = PoM^ , (45) 
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where f refers to the Hermitian conjugate. This results in the detailed balance 
for transition probabilities, and vice versa, when those are introduced. 

By using the constraint Eq. (41) one has 

N N 

MQVo^ -kBllVZ^QVo ^-kBY.Y.'^n® i?^"'^^ = (46) 

n=0 p=0 

which proves that the equilibrium density superoperator is a right eigenfunc- 
tion of the master supermatrix with an eigenvalue zero. Here, we have used 
the fact that 

So that, Eq. (43) can be rewritten as 

^^-MQV. (48) 

Now, by substitution of the master equation (48) into the generalized Liouville 
equation (20) we obtain 

^V{t) + ih-^VC V{t) = -QM V{t) , (49) 
where QM. should have the same structure as QC, thus 

N 



QM= J2 P^in+D^n ® Tr { M("+i|^') } (50) 

and 



n,p=0 ^ ' 



N 



QMQV=J2 Wi)^n® Tr {m("+i|p)0Z}(p)} . (51) 

n,p=0 ^"+^) 

This last equation allows us to write 



+ ir' {hSd%> - H^;^D^) = (52) 

equivalent to Eq. (21). 

Here, Eq. (52) can be put into a more common form by introducing the tran- 
sition superoperator wj^^^^ defined through 

•^^^ a/3,75 — ^^^i5,al3 + -'*-a/3,7<5 ' K"^^) 

where 

Aiti = Ia,aM;^Kp . (54) 
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The auxiliary matrix "^^g is not arbitrary, satisfying 

EW^iS = <^7^*i7^<^nP> (55) 
a 

rin\p) 

Y detailed balance, A^*^^^^ should t 

one has 



in view of Eq. (38). This means that '^aa% ^ diagonal matrix. In order 

i\p) 



to satisfy detailed balance, A^"^^^^ should be Hermitian. Hence, with Eq. (45) 



constituting the expression of the detailed balance principle. Thus, projecting 
Eq. (52) onto its diagonal part by means of the tetradic operator Pij^ki — 
SijSik5ji [19], 

P (^^^^"^ + [h^^'K ^^"^] j = (57) 

and by using Eqs. (53)-(55) we obtain 

+ ^f^-' (hSd^^J - H^^D^-J) = (58) 



Qf7 7Q a7 7tt 



For n — 1 

^D^l + K7 - ^S^S) = (59) 
Trlz^^'.V^f) _ fy- 1^(212) ^ ^(2)1 _ 

1^ 70 70, aa ^ / ^ 77, aa y aa j 

By neglecting the correlations which means assuming D^^j (12) = D^^j (l)L'^^j (2) , 
we obtain 

|^«(i) + ^^-^ K^^Sd) - <^<ni)) 

- Tr {/^iy(l)i^ii^(2)<ti(12) - (e W^S.(12)) i^i^i(l)i^«(2)} (60) 
which constitutes the quantum Boltzmann equation. 
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5 Conclusions 



Here we have shown a representation of the statistical description of an N-body 
system in terms of a density superoperator which reveals the irreversibility at 
the macroscopic level. In the scenario drawn by the Liouville-von Neumann 
equation with its associated Gibbs-von Neumann relative entropy, irreversibil- 
ity is hidden because of the inherent unitary evolution of the N-particle density 
operator. Nonetheless, in the framework of the quantum BBGKY-hierarchy 
the density superoperator follows a non-unitary dynamics given by the gener- 
alized Liouvillc-von Neumann equation. Thus, the generalization of the Gibbs- 
von Neumann entropy in terms of the density supoeroperator we propose is 
not conserved by this dynamics. Moreover, this dynamics is trace preserving 
which as is known [14]- [18], ensures the non- negative character of the entropy 
production which is coherent with the second law of Thermodynamic. 

In addition, following the procedure of the non-Equilibrium Thermodynamics 
we introduce a master matrix which satisfies a kind of Green-Kubo relation. 
This, allows us to transform the BBGKY-hierarchy into a hierarchy of master 
equations which we rewrite in terms of a transition superoperator related to 
the master matrix. Once this has been done, by neglecting correlations in the 
spirit of the Boltzmann's Stossahlansatz, we derive the quantum Boltzmann 
equation. 

To summarize, our point is that the framework of the BBGKY-hierarchy which 
explicitly incorporates the correlations between all the particle clusters in the 
system constitutes the most appropriate description leading to irreversibility 
at the macroscopic level. To manifest this irreversibihty all that is needed is 
a coherent definition of the entropy which is our generalaized Gibbs-von Neu- 
mann relative entropy in addition to the generalized Liouville-von Neumann 
equation. 

In this way a macroscopic time arrow arises in a natural way in the system. 
Hence, we can say that our analysis constitutes a microscopic basis of the 
macroscopic irreversibility. 
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